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Abstract 

Based on the "cannonball model" for gamma-ray bursts of Dar and De Riijula it is 
proposed that masses of baryonic plasma ( "cannonballs" ) , ejected in bipolar super- 
nova explosions in our Galaxy are the sources of hadronic Galactic cosmic rays (CRs) 
at all energies. The propagation of the cannonballs in the Galactic disk and halo is 
studied. Two mechanisms for the acceleration of the observed CRs are proposed. The 
first is based on ultrarelativistic shocks in the interstellar medium and could accelerate 
the bulk of CRs up to the "knee" energy of 4 x 10 15 eV. The second operates with 
second-order Fermi acceleration within the cannonball. If the total initial energy of 
the ejected plasmoids in a SN explosion is 10 53 ergs or higher this second mechanism 
may explain the CR spectrum above the knee up to the highest observed energies. It is 
shown that together with plausible assumptions about CR propagation in the Galactic 
confinement volume the observed spectral indices of the CR spectrum can be theoret- 
ically understood to first order. The model allows a natural understanding of various 
basic CR observations like the absence of the Greisen-Zatsepin cutoff, the anisotropy 
of arrival directions as function of energy and the small Galactocentric gradient of the 
CR density. 



1 Introduction 

1.1 Motivation for a "unified" (single-source class) explanation of the 
spectrum of Galactic cosmic rays 

The flux of the dominating hadronic component of the local non-solar cosmic rays (CR), 
has been measured between about 0.1 GeV and 3 • 10 20 eV. Neglecting solar modulation 
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effects, its energy spectrum can be well described by a single power-law that changes its 
slope slightly at only two energies - at ~ 4 • 10 15 eV (the "knee") and ~ 3 • 10 18 eV (the 
"ankle" )[]. This striking simplicity and unity of the data over more than 10 decades in 
energy originally led most authors to ascribe the origin of hadronic cosmic rays to a single 
"universal" source class. 



In an influential review, Ginzburg[31| suggested that the remnants of supernova explosions 
in our Galaxy (SNRs) are the dominant source of CRs at all energies. However, SNRs 
cannot explain the origin of CRs with energies above above about 10 18 eV because at these 
energies charged particles are certainly no longer confinable in the relatively weak magnetic 
fields of SNRs. In another classical paper Burbidge ||i~0| proposed the radical alternative 
of an extragalactical origin of CRs at all energies. In this scenario it proved difficult to 
understand the origin of CRs with low energies: no extragalactic source class seems able to 
fill up the universe with the rather high total energy-density of the locally observed CRs 
(~ 0.5 eV/cm 3 ) within a Hubble time. 

As a compromise, an "eclectic scenario" for CR origin was proposed by Morrison pS|]. This 
idea became widely accepted in the 1970s in the following form: SNRs accelerate the CRs 
below the knee and extragalactic sources are the source of CRs above the ankle. Between 
the knee and ankle another class of objects - of an as yet unclear nature - is supposed 
to accelerate the CRsQ. However, the ankle lies just at the energy where CRs become 
unconfined from the Galaxy. This fact can only be causally understood if the same universal 
source class supplies the cosmic rays below and above the ankle [66, |53[] . The eclectic scenario 
was born of distress, namely of the difficulty that no "universal" source class was known 
that can accelerate particles to the highest observed energies and at the same time supply 
the rather high local energy density of CRs. The aim of this paper is to revisit the idea 
of a universal origin of Galactic CRs and to propose a source class capable of accelerating 
nuclei to the observed energy spectrum. 



1.2 Clues to the nature of the "single source" class from experimental 
evidence on cosmic rays 

In all probability the universal source class resides within or near to our Galaxy. It is all 
but impossible to supply the high ambient Galactic cosmic-ray energy density with any 
extragalactic source class. The absence of the Greisen cutoff in the experimentally deter- 
mined cosmic-ray spectrum is another strong argument in favour of a Galactic source class. 
The observed lack of an overall anisotropy of arrival directions in UHE cosmic rays, the 
very small anisotropy at all smaller energies, and the - at best - very weak concentration 
of the cosmic-ray density towards the Galactic centre forces us to conclude that the "uni- 
versal source class" (if it exists) has a very different spatial distribution than SNRs which 

1 Possible further small deviations are discussed in section 2.2.1 

2 Biermann's[B] theory of CR origin does not make this assumption and is further discussed in section [L4 
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Figure 1: A scheme of the scenario proposed here. Energetic core-collapse supernovae in 
the Galactic disk and centre eject plasmoids ( "cannonballs" ) that are slowed down to non- 
relativistic speeds in the Galactic halo. These cannonballs are proposed to accelerate the 
major part of hadronic CRs at all energies. 



populate the disk and central region of the Galaxy. It seems to reside in the Galactic halo. 
A possible origin of CRs at all energies in the Galactic halo, due to Fermi acceleration in 
hydromagnetic turbulence, has been discussed by BurbidgeQ. Biermann & DavisQ and 
Plaga[53] discuss mechanisms by which a large Galactic halo may be filled homogeneously 
with CRs. 



1.3 Ejecta from bipolar supernovae as sources of Galactic cosmic rays 

There has been recently strong observational evidence that core-collapse supernovae (SNe) 
can produce "long" 7-ray bursts (GRBs) ^(|. A related recent revolution in the understand- 
ing of SNe is observational evidence that apparently all core-collapse SNe are asymmetric 
"bipolar" explosions ]69[. Observations cannot yet directly decide yet if: 

a. only very few "exceptional" SNe emit GRBs[|5l|] or 

b. a sizeable fraction of or even all core-collapse SNe emit GRBs jet like into a narrow cone 
along a symmetry axis\ 11]. 

There is direct observational evidence that GRBs are emitted by matter at ultrarelativistic 
speeds. Circumstantial evidence suggests that GRBs involve rather strong magnetic fields. 
This suggests that they are excellent candidates for the acceleration of CRs at all energies. 
Especially if option b. above should be correct, the following question will be decisive for 
CR physics: 
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Do the well known non-relativistic isotropic ejecta (SNRs) or the only recently 
discovered relativistic jet-like ejecta dominate the CR production of all super- 



novae : 



It has been suggested that the very high ejection speeds in conjunction with strong magnetic 



fields allow acceleration to much higher energies than SNRs [68, 7C]. This makes jet-like SN 
ejecta serious candidates for the long sought universal source class of hadronic CRs at all 
energies |L7|, 

The potential importance of Gamma-Ray bursters for the origin of Galactic CRs at ener- 



gies above the knee had been first studied by Milgrom and Usov[48] and Dar|l?]. Recently 



Pugliese et al.[56] concluded that GRBs are unimportant as sources of CRs at all energies 



whereas Dermer[22] argues that objects related to GRBs might supply most of the observed 
low-energy hadronic CRs. The reason for this disagreement is the use of very different GRB 
models. 



1.4 Aims and plan of this paper 

The basic idea of the present paper - that ultrarelativistic ejecta from stellar-collapse events 
accelerate the dominant part of hadronic CRs at all energies and release it in the Galactic 
halo (fig. |T|) - was first proposed by Dar & Plaga]l9[|. At that time the nature of the GRB 
emitting collapse-events was unclear (supernovae were only mentioned as one possibility) 
and they were called "Galactic Gamma-Ray Bursters (GGRBs)". In the present paper the 
nature of these collapse events, the propagation of their ejecta and the CR-acceleration 
mechanisms operating near them are discussed more quantitatively. 

Section § discusses what one can conclude about the properties of a hypothetical "single- 
source class" purely from the phenomenology of CRs. I conclude that universal sources 
must be mainly located outside the solar circle, in the Galactic halo. 

Further observational evidence in favour of such a location of the CR sources is discussed 
in section [|. The Galacto-centric distribution of 7-rays and anisotropy of charged CRs 
support the idea that CRs below the ankle have their origin mainly (but not exclusively) 
at Galacto-centric distances exceeding the solar one. 

A candidate for the "single-source class" is proposed in section ^ ultrarelativistic plasmoids 
ejected in bipolar supernova explosions ( "cannonballs" (2(|). The motion of cannonballs 
in the Galactic disk and halo is discussed. Section || studies mechanisms with which these 
ejecta produce the hadronic CR spectrum at all energies (||). 

This scenario is related to ideas of Biermanj?]] who proposes an origin of Galactic CRs up 
to the ankle in SNRs and an origin of the higher energy CRs in extragalactic jets. My work 
can be seen as a proposal to "replace" the non-relativistic SN remnant with ultrarelativistic 
supernova ejecta. These ejecta ("cannonballs") "assume" the role of an extragalactic jet in 
a later stage of their evolution. 

The ideas of sections [2] and ||| are logically independent. It is conceivable that there is 
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a single-source class, but I err in its identification, or that cannonballs contribute to the 
CR spectrum only in certain limited energy ranges. Section ^ confronts expectations in a 
scenario where cannonballs dominate the CR production at all energies with observations. 
Section |7] summarises the main assumptions, achievements and predictions of this paper. 
This paper only treats hadronic CRs, leaving aside the origin of electrons in the Galactic 
CRs as a completely separate issue. 



2 Properties of a hypothetical single-source class of cosmic 
rays derived from observations 

In this section I try to derive the likely properties of a hypothetical "universal" CR-source 
class from observed CR properties without making assumptions about the physical nature 
of the sources. In this I optimistically take it for granted that a "simple" understanding of 
the CR spectrum - e.g. without accidental cancellations of additional effects - is possible. 

2.1 Energy dependence of cosmic-ray diffusion coefficient assumed in this 
work 

In this work I will assume that the cosmic-ray propagation is purely diffusive. The depen- 
dence of the diffusion coefficient D on the total energy per nucleus E is usually parametrised 
as: 

(1) D ~ E a 

There are two theoretically motivated values for a. A plasma with turbulent inhomo- 
geneities following a "Kolmogoroff spectrum" would lead to a diffusive CR motion with 
a=0.33, alternatively, with a " Kraichnan spectrum" a=0.5 would apply[12|. Direct de- 



terminations of the electron-density fluctuations in the interstellar medium indicate a 
turbulence spectrum that leads to a power-law dependence of D according to eq. (JT[) up to 
very high energies^]. However, such measurements are not precise enough to determine a. 
CR data obtained at low energies indicate a m 0.6-0. 7j5|] . This value can be made com- 
patible with both theoretical values by postulating varying degrees of "re-acceleration" 
during the propagation of cosmic rays. There is recent experimental evidence for some 
re-acceleration [16] but it is not yet clear which value for a is preferred by these data. 



It is generally accepted and plausible that CRs are magnetically confined to the Galaxy up 
to energies of about E ank i c ps 3 x 10 18 eV|]. At lower energies the power-law index of the 
CR spectrum changes only once by 0.3 at the knee. In principle the "knee" could signal a 
change from a=0.33 to a=0.5. However, this is unlikely because the observed change in the 
power-law index of ~ 0.3 is significantly larger than the difference between the theoretically 
preferred indices of (0.5 - 0.33). I conclude that in any simple scenario - where a has one 

3 up to scales corresponding to the Larmor radius of ~ 10 17 eV particles in the Galactic magnetic field. 
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of the theoretically preferred values at all energies - a is constant at either 0.33 or 0.5 at 
all energies up to the ankle. In the next section I will argue that - at least in a scenario 
with "universal sources" - a=0.5 is the preferred value. 

2.1.1 An argument in favour of a=0.5 at all energies valid under the assump- 
tion that universal sources of CRs exist (can be omitted on first reading) 

Let d be the typical distance of the universal sources from the edge of the Galactic confine- 
ment volume (see fig.(|l])). The CR flux at an energy of 

Ei = 40 GeV is a factor K=(Ei/E an kj e ) _Q higher than at energy E an kj e due to a K times 
longer confinement time r. The confinement time at energy Ei is given as t\ ~ d 2 /Di, 
where Di is the diffusion coefficient at energy E^. For a halo size h ~ 20 kpc (a large halo 



size is preferred for for the universal-source scenario, see section 2.3.2) D was determined 
from measured isotope-ratios as approximately 9 x 10 28 cm 2 /sec[33|. The confinement 
time at energy E an kj e is approximately r a « d/c, the value for unconfined motion. Setting 
K=ri/r a I obtain: 

d « (^l/^kle) -Q X ZVc 

~ 8 kpc for a = 0.5 
(2) ss 0.3 kpc for a = 0.33 

The confinement volume is certainly extended by more than 10 kpc in our scenario (see 
section 2.3.2j ) and therefore the case a=0.33 would lead to a "crowding" of the universal 



sources in the outer < 3 % of the confinement volume. This seems rather unnatural and 
we therefore assume in the following that a = 0.5 at all energies < E a . 
To summarise, the assumption of a E 0,5 dependence of CR diffusion at all energies where 
Galactic confinement occurs is theoretically motivated, in agreement with all data and 
allows a natural spatial distribution of universal sources. 

2.2 The value of the spectral indices of the source spectrum, before 



Galactic modulation( [|T9|| ) 

The assumption of a=0.5 (eq.||) fixes the spectrum of the "universal CR sources". Let the 
source spectrum be given as: 



(3) ^source ~ E S 

If this source emits into a volume that diffusively confines CRs the observed spectrum 
within this volume is given as: 

(4) F obscrvcd (£ < £ anklc ) ~ ~ 

4 The diffusion coefficients are specified by Strong & Moskaienko (1998) at a rigidity of 5 GV. For 
intermediate masses (A=16) this corresponds to an energy per nucieus Ei « 40 GeV. 
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This spectrum holds at all energies where cosmic rays are magnetically confined to the 
Galaxy (i.e. up to the ankle). At higher energies the pure source spectrum is expected 
again: 

(5) -^observed (E > -Eanklc) ~ E 5 ~ E 13 

at higher energies. The observed index 72 between knee and ankle is {82 + a)=— 3. I infer 
for the source index 82=— 2.5, a value that is compatible within 1.3 a with world- mean of 
experimental determinations of the index above the ankle of 73=— 2.75 ± 0.2|50|] as it must 
be within my model. According to section |2.1| I assume no breaks in the energy dependence 
of the diffusion coefficient which governs CR propagation, so the "knee" in the CR spectrum 
must be a source feature. The spectral index of the source must then decrease by 0.3 at 
the knee, so the 8\=—2.2 at energies below the knee. This general scenario for the origin 
of the total CR spectrum is illustrated in fig. (|3|) (taken from Ref . 0] ) . 



o 
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Figure 2: A schematic sketch of the main steps in the scenario discussed here. The shaded 
balls symbolise the cannonballs. The thin un-arrowed line in the second image stands for 
a shock wave driven by the cannonball into the interstellar medium. Step 1 is further 



explained in section 4.1, step 2 in section 5.1 and step 3 in section p 



2.2.1 Unsettled issues: structure at the "knee" and the "second knee" 

The discussion of the spectrum in the previous section is oversimplified. Erlykin and 



Wolfendale|2£] have argued that an additional component that raises the flux by about 
a factor 2 in a narrow energy interval near the knee at 4 x 10 15 eV is indicated by the data. 
Detailed models for the "universal sources" have to deal with this feature. 
At energies around 4 x 10 17 eV the spectrum steepens by about 0.2 in the power law in- 
dex, thus forming a "second knee" |)TJ] . In our scenario it seems natural that the transition 
between diffusion in a turbulent plasma and approximate straight-line motion is gradual 
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and therefore a steepening occurs before the ankle. 

I defer more detailed speculations about the nature of these finer details of the spectrum 
to a later time. 



2.3 Spatial location of the universal CR sources 

2.3.1 Upper limit on distance of universal sources from Galactic centre: iso- 
topic data of low-energy CRs 

The data on isotopic composition of CRs yield an upper limit on the size of the Galactic 
halo in which magnetic confinement takes place. Because the CR spectrum can only arise 



in the way outlined in section 2.2 if the particles are initially magnetically confined, the 
maximal size of the halo is also an upper limit on the distance of the universal CR sources 
from the Galactic centre. 

The size h of a cosmic-ray halo around the Galaxy is only weakly constrained by results on 
the abundance of unstable isotopes in low-energy cosmic rays fo3|j . A halo with h > 4 kpc 
is required to exist, and h 3> 20 kpc seems to be ruled out with all parameter choices, so 
we assume a maximal allowed size of 30 kpc^]. 

2.3.2 Lower limit on distance of universal sources from Galactic centre: isotropy 
of UHE CRs 

No significant anisotropics have been found in the sky distribution of CRs with energies 
above 4 x 10 19 eV[5C, [73|]. However, with a world statistics of 114 events, anisotropies on 



the order of up to 30 % amplitude cannot be ruled out. At these energies - far above the 
"ankle" - light nuclei are expected to move in straight lines to good approximation in our 
scenario. Also light nuclei comprise a major fraction of all CRs in "universal scenario" . (fig. 
|8|). We conclude that a major part (> 50 %) of all universal sources must lie at distances 
beyond the solar circle (i.e. Galacto-centric distances > 10 kpc). 



2.3.3 The likely location of "universal" CR sources 

The above consideration indicate a typical location of the "universal sources" at about 10 

- 30 kpc from the Galactic centre. For the reasons mentioned above these limits (especially 

the upper one) are quite uncertain. A detailed simulation using the code of Strong & 

5 This upper limit is quite uncertain mainly for four reasons. 1. Strong and Moskalenko Js3[] do not 
take into account uncertainties of the input parameters, which are considerable. 2. The spatial distribution 
of universal CR sources is expected to be different from the one assumed by Strong and Moskalenko. 3. 
CR-electron sources might have a very different spatial distribution (following the SNR distribution) from 
the one of CR-hadron sources (following the cannonball distribution in our scenario). 4. The plausible 
possibility of a diffusion coefficient that is smaller in the Galactic disk than the haloj32| is not considered. 
This possibility gains in potential importance if sources are not located in the Galactic disk. 
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Moskalenko|63|] and taking into account the special features of the universal source class is 
urgently required. 




Figure 3: The observed flux of cosmic rays (thin line) as a function of primary energy E is 
well described by a power law that changes its slope sharply at only two energies, the "knee" 
and the "ankle". At energies below the ankle it is enhanced (by a factor (E/E an fc/ e )~ ' 5 ) 
over the Galactic Gamma-Ray Burster (GGRB) source spectrum (thick line, a power law 
with differential power law index of -2.2 below the knee and ~ -2.5 above it) by way of 
trapping in the Galactic halo magnetic fields. The "GGRBs" are identified with bipolar 
SNe in the present paper. 



3 Observational evidence in favour of CR sources located 
mainly (but nor completely) outside the solar circle 

3.1 Galacto-centric gradient of cosmic rays from 7-ray data 

The analysis of EGRET data from Mattox and Strong[62] and Erlykin et al.[^5| find that 
the cosmic-ray density is constant within about 25 % between 5 and 15 kpc (5 - 18 kpc in 
Erlykin et al.(25|]) distance from the Galactic centre. In this range the SNR- and pulsar- 
densities density decrease by about a factor 5 to 10. Were these object classes the main 
sources of CRs - and taking into account CR diffusion in a large halo - a much larger 
decrease of CR density by a factor > 3 over this range of Galacto-centric distances would 
be expected |8]]. The last bin (15 to 30 kpc in Mattox and Strong^]) shows a density about 
30 % lower than the previous ones. The analysis of the EGRET team [[37]] shows a drop in 
the density by about 30 % in the range 5-10 kpc than a rise in density up to 15 kpc in 3 
quadrants. At larger Galacto-centric distances a precipitous drop in density was inferred. 
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These results are not in complete agreement with each other but they do suggest that the 
universal sources of Galactic CRs are mainly located in shell at about 20 kpc distance 
from the Galactic centre. In this approximation one expects a constant density at smaller 
distances and a rapid drop at larger distances; this is in agreement with the data. 

3.2 Anisotropy of arrival directions of Galactic cosmic rays 

If all universal CR sources were located outside the solar circle one would expect no 
anisotropy of Galactic cosmic rays at any energy: there is none inside a shell isotropi- 
cally emitting cosmic rays. There is however experimental evidence for small anisotropies 
in Galactic cosmic rays at all energies p^] . The observationally well established energy- 
independent anisotropy with an amplitude of ~ 8 x 10 -4 at energies between about 1 and 
100 TeV is probably the result of a motion of the local rest frame of cosmic rays relative 
to the solar system with a speed of about 80 km/sec via the Compton-Getting effect |l5|] . 
The nature of this motion is unclear. SchmelejpO]] recently proposed an anomalous velocity 
component of the solar-system surroundings connected with a motion of the "local fluff" 
out of the centre of the local super-bubble in which the CRs are at rest. If this interesting 
suggestion is correct, the anisotropy below 10 14 eV contains information about the plasma 
that confines CRs but holds no direct clues about the origin of cosmic rays. 
At higher energies there is a large world-data set which indicates an anisotropy that in- 
creases with energy E as E 0,5 up to the highest energies (fig. ||) |2JJ]. In spite of a very large 
systematic scatter, the measurements of the direction of first-harmonic maximum intensity 
in the total world-data set seem to cluster around the Galactic-centre direction (see fig.||) . 
Because individual anisotropy measurements usually have marginal statistical significances 
and the derived locations of maximum intensity do not agree well, doubts have been ex- 
pressed if these results are more than upper limits on a possible anisotropy. However, a 
recent measurements and reanalysis indicate significances for anisotropies above 5 a. 
Clay et al.jl4|] reanalysed the world data set in the southern hemisphere in the energy range 
10 15 - 3 • 10 15 eV and found an anisotropy amplitude of 

(6) S a = 0.33 ± 0.06% 

which corresponds to a 5.5 a effect. The direction of maximum intensity has a right as- 
cension of 19.7 ± 0.7 h, not far from the one of the Galactic centre (at 17.8 h). Hayashida 
et al.|3^] report an anisotropy at energies above 10 18 eV of 3.9 ± 1% with a statistical 
significance of 5.4 a. Again an excess mainly from the Galactic-centre direction is indi- 
cated. These two high-significance results are in reasonable agreement with each other if 
the E 5 energy dependence of the anisotropy amplitude of the world data set is assumed. A 
more detailed future analysis needs to take into account the different geographical locations 
of the various experiments which necessitate small corrections of the measured anisotropy 
amplitudes and directions. 
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Figure 4: World-data set as of 1997 of measurements the anisotropy amplitude of the first 



harmonic as a function of energy. From the thesis of Schmele[60] where all references are 
given. The thin dotted lines indicate the data point obtained in this thesis from data 
obtained by the HEGRA air-shower array. The diamonds indicate the recent results by 
Clay et al.[|l4|] (near 1 PeV) and Hayashida et al.|33|] (near 1 EeV), discussed in the text. 
The quoted errors of these two determinations are smaller than the symbol sizes. The thick, 
full line indicate the anisotropy theoretically predicted (eq.(|?]) with /=0.15 (section |6 . 2 . 1| ) ) . 
The thick dotted line indicates the constant anisotropy at energies between about 1-100 
TeV probably due to an anomalous motion of the solar system relative to the rest frame of 
CRs (see text). 



Summarising, two new high significance measurements confirm the long known indication 
of an anisotropy amplitude rising with square root of energy and prefer an origin of this 
anisotropy of arrival directions from the direction of the Galactic centre. 
A natural origin for the observed anisotropy is a net outward-streaming motion of cos- 
mic rays produced in the Galactic centre region. If such a streaming motion exists, the 



Compton-Getting effect [15] leads to a large scale anisotropy with a maximal intensity near 
the Galactic centre. Let us estimate the expected magnitude an anisotropy 5 a produced in 
this way. The preferred values of the energy dependence of the diffusion constant D, halos 
size h and spectral index 7 as discussed in sections |2.1.1| and 2.3 were chosen. Here "/" is 
the fraction of CR sources that are located within the solar circle. 

5 a = 0.5%(£>/9 • W 28 cm 2 /sec)/(h/20kpc)/c x 
(7) ((t/2.7) + 1)(//0.1) • {(E/10 15 eV)/A0GeVf 5 

This expression is valid from about 10 15 eV (at lower energies another source of anisotropy 
dominates) to about 3 x 10 18 eV (where CRs become unconfined to the Galaxy). Here "/" 
is the fraction of CR sources that are located within the solar circle. Eq. (|7|) is plotted in 
fig.0 and is seen to agree with the data within the rather large errors. 
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Figure 5: World-data set as of 1997 of measurements the direction of the maximum of the 
CR intensity in right ascension as a function of energy (first harmonic) . From the thesis of 
Schmele[|6C|] where all references are given. The thin dotted lines indicate the data point 
obtained in this thesis from data obtained from the HEGRA air-shower array. The thick, full 
line indicates the direction towards the Galactic centre. The diamonds indicate the recent 



results by Clay et al.|14] (near 1 PeV) and Hayashida et al.|3^] (near 1 EeV) discussed in 
the text. The error bars of the result of Clay et al. (1998) are comparable to the symbol 
size, no error was given in Hayashida et al.|3^] for the determination of the direction. The 
thick dotted line indicates the range where an energy-independent anisotropy was found 
(see text). 



A consistent explanation of the CR anisotropies at all energies is thus possible 10 % 

of the universal sources are located in the Galactic centre and 90 % outside the solar circle. 



This is consistent with the observational constraints of section |2.3| and the prediction within 
the cannonball model (section |6.2.1| ). 



4 A candidate for the "universal" CR sources: ultrarelativis- 
tic ejecta from bipolar supernovae (aka "cannonballs" ) 

4.1 Ultrarelativistic ejecta ("cannonballs") from supernovae 
4.1.1 The cannonball model of Gamma- Ray Bursts 



Recently Dar and De Rujula[20, |2l[ developed the "cannonball model" for Gamma-Ray 



Bursts (GRBs). Similar to Woosley's[75] "failed supernova" concept it proposes that the 
mechanism for core-collapse supernovae can involve a phase where matter falls onto a cen- 
tral compact object from large distances. Episodic accretion onto this central black hole 
or neutron star leads to the subsequent repeated emission of distinct masses of plasma, the 
"cannonballs" (fig.0), initially moving with a Lorentz V ~ 1000. Radiation emitted by the 
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cannonballs is quantitatively shown to explain the properties of GRBs. This radiation is 
emitted isotropically in the rest frame of the cannonball, but in the observer frame the 
Lorentz boost strongly collimates the radiation into a cone with opening angle 9 =l/r 
along the direction of motion ( "relativistic beaming" ) . Therefore only a very small fraction 
~ # 2 /4 ~ 2 x 10~ 7 of all cannonballs are visible. 



Recently I found [ 54 that it is possible to quantitatively understand a "cepheid-like" re- 



lationship between variability and total luminosity of GRBs [27] and between the width 
and height of spikes within a given GRB|58| within the cannonball model. This makes it 
plausible that the broad features of the cannonball model are correct. I will assume that 
this is the case in the rest of this paper. 

The initial energy of these cannonballs is not well known. If the core of a pre-collapse star is 
rapidly rotating, 3 • 10 53 erg may be stored as rotational energy of the collapsed object, 



and the jet seems the most effective way to dispose off this energy 1 52]. I will assume below 
a total energy of 10 53 erg ejected in each of the symmetric cannonball jets in a SN (rather 
higher than what is normally assumed but in line with the expectations of Dar and De 
Rujula pl|| ). Smaller values lead to conditions in which CRs cannot be accelerated to the 
observed ultra-high energies (section 



4.1.2 The frequency of supernovae that eject cannonballs, inferred from the 
observed GRB rate 

For some GRBs Fenimore and Ramirez-Ruiz [27] infer redshifts in excess of 12 from their 
variability - luminosity relation. Because it seems unplausible that much larger redshifts 
occur, within the cannonball model this indicates that cannonballs ejected exactly into 
the observer direction are visible from the entire universe. If that is true, the rate of 
cannonball-ejecting supernovae in our Galaxy is given as [|D| : 

50years lCr/year 

(-ill v^mea.n 

17 ^(r/iooo) 2 )' 

((L B (Galaxy)/2.3 x 1O 1O L )/ 
( PB /1.8hl0 8 LoMpc- 3 ))/ 
x 15 



(8) (#sfr(z=o)/( /(I + z)- 1 R SFR(z) (dV c /dz)dz)) 



Here -Rgrb is the observed total rate of GRBspq], Lb is the B-luminosity of our Galaxy, 
Pb is the B-luminosity density in the local universe |47j, Rsfr is the star formation rate 



at a respective redshift derived from optical observations [42]. The preferred value was 
derived with a volume element dV c for Q=l and A=0. #mean * s the mean squared angle 
between an observed GRB and the observer divided by T 2 . This angle is expected to be 
somewhat larger than 1/T 2 because nearby GRBs are bright enough for detection also at 
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viewing angles > 1/T. The mean(# 2 ) ~ 5 was inferred by calculating the O 2 for all bursts 



in Fenimore and Ramirez- Ruiz p7| using relations valid in the cannonball model[54]. An 
isotropic luminosity at 9 = of 2 x 10 55 erg/sec (the highest luminosity of any GRB in 
their sample) was assumed for this estimate. The preferred value of -R pr0 g = 1/50 years 
corresponds to the estimated total rate of all core-collapse supernova in our Galaxy f59fl. 
thus indicating that all core-collapse supernovae eject cannonballs. 

The major uncertainty in the estimate of R pr0 g is the dependence on the star-formation 
rate (SFR) on z. In particular Fenimore and Ramirez- Ruiz p7| derive a -Rsfr-grb based 
on observations of GRB space densities alone. This SFR continues to rise at redshifts 
beyond 1, in contrast to results based on optical observations. Using these results I obtain 
^SFR-GRB( z =o)//( 1 + z)' 1 RsFR-GRB(z)(dV c /dz)dz « 250. This corresponds to R pmg « 
1/800 years. The possible choice of other world-models introduces another uncertainty of 
about a factor 2 in this estimate, so the plausible range of rates of cannonball-producing SNe 
is about 1/50 - 1/1600 years. Preferred candidates for jet-forming core-collapse events have 
similar formation rates: SN Ib/c supernovae (R ~ 1/300 years p9[|), "failed supernovae" (R 
» 1/1000 yearsf|] ) and "collapsars" (R » 1/1500 yearsf29j). 

Summarising, based on the observed GRB rate, the progenitor rate for cannonball ejecting 
supernova could lie between about 0.05 - 1 of total core-collapse SN rate. 



4.1.3 The total cannonball-energy converted into CRs 

The "CR luminosity" L cr - the energy per time in form of hadronic CRs injected into the 
Galactic confinement volume - is given as: 

L cr = 2 x 10 42 erg/sec(-Rp rO g/1000years) x 
(9) CEW(2 x 10 53 erg))(e/0.33) 

Here e is the conversion efficiency of total energy to hadronic CRs within the collimated 
(section |4.2j ) cannonball, that is likely near the equipartition value of 1/3 in our scenario 
(see section |5.2| ). i? pr0 g is the rate of cannonball-ejecting supernovae in our Galaxy. If 
cannonballs are the dominating source of hadronic CRs they must supply the observed 
"CR luminosity" L cr (exp) ~ 10 41 erg/sec which can be inferred from experimental data|J. 
Within the plausible range for i? pr0 g discussed in the previous section, less than 1/40 of 
the total cannonball-energy has to be converted the energy of CR particles to explain the 
observed CR luminosity L cr (eq.(||)). 

4.2 Collimation of the cannonball 

4.2.1 Assumption of cannonball confinement 

In the model of Dar and De Rujula|2(J the cannonball initially expands with the speed of 
sound in a relativistic gas c/y/3. It is a basic hypothesis of the present paper that some colli- 
mation mechnism prevents any further expansion of the cannonball once the cannonball has 
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reached a limiting radius iiWiit ~ 0.3 pc. If there is equipartition between magnetic-field 
energy density, turbulent plasma motion and cosmic-ray energy density a total magnetic 
energy of E tot /3 ~ 3 xlO 52 erg (see sect. |4.1.1| ) corresponds to a magnetic field inside the 
cannonball of B = ^ E tot /Rf imit & 0.6 G. 

Efficient collimation is commonly observed in extragalactic jets (for a review see pSf ). There 
is no complete theoretical understanding of the mechanisms leading to collimation in gen- 



eral. Longair recently wrote [45]: "By which means relativistic jets are collimated, is one of 
the major unsolved problems of astrophsics. Almost certainly, the mechanism will require 
the presence of magnetic fields...". In principle a slight deviation from equipartition with 
strong magnetic fields might serve to contain turbulent plasma and CRs. However, the 
"plasma-virial theorem" forbids a purely magnetic stress balance without additional outer 
stresses^]. The only viable candidate for the outer stress in the case of cannonballs seems 
to be ram pressure from the ambient interhalo medium. It has been shown long ago that ram 
pressure can confine a spherical mass of plasma f23|] , and Heinz recently argued that the same 
mechanism might lead to confinement of "bullets" ejected by GRBsp5|. Equating the typ- 
ical expected ram pressure onto a subrelativistic CB[] P ram = rn p (ni n t e rhaio/(10 -3 /cm~ 3 )) 
(r/1.2) c 2 with its internal magnetic-field pressure B 2 nt /47r shows that with the parameters 
indicated above ram pressure can only balance an internal field Bj n t of about 10 mG, about 
a factor 50 smaller than required. Efficient collimation therefore requires a magnetic field 
that increases towards the interior of the CB (similar to what happens in a tokamak). What 
magnetic-field configuration of the cannonball would allow a stable confined configuration 
in conjunction with external ram pressure is an open question, closely related to Longair's 
"open problem" mentioned above. As a detailed simulation of the internal dynamics is 
beyond the scope of the manuscript, below I make the simplest possible assumption of free 
expansion until a limiting radius Rii m it is reached. When the cannonball has slowed down 
to nonrelativistic speeds and the ram-pressure becomes very small, a final Sedov- Taylor 
expansion must take place pC|], 

4.3 Motion of confined plasmoids in the Galaxy 
4.3.1 The interstellar medium ambient to cannonballs 

Lingenfelter, Higdon and Ramaty] ^] have pointed out that about 85 % of all core-collapse 
supernovae are expected to occur in "super-bubbles" . These are cavities in the interstellar 
medium (ISM) blown by repeated SN explosions and filled with hot (10 6 K) and tenu- 
ous (number density 10 _3 /cm 3 ) gas. Their conclusion seems compelling to me. It is an 
undisputed observational fact that about 50 % of the ISM volume of the Galactic disk is 
filled with super bubbles |H|]. Therefore at least half of all supernovae are expected to 
occur in super bubbles. The actual fraction must be somewhat higher than this because 



see discussion in section 5.2 of the CB propagation 
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core-collapse supernovae are known to be spatially and temporally correlated. 
The typical super-bubble size is several hundred parsecs, and a large fraction of the super- 
bubble borders on the Galactic-halo medium. The Galactic halo is observationally known 
to contain hot, tenuous plasma which has similar properties (and possibly the same origin) 
as super-bubble medium. At a distance of 50 kpc from the Galactic centre the halo density 



was roughly determined as 10 /cm [71] and from absorption studies the halos of spiral 
Galaxies are known to extend to several hundred kpc] 41 1. 

To good approximation the cannonballs are thus expected to propagate in 85 % of all 
cases in super-bubble and halo medium and in 15 % of all cases in the normal interstellar 
medium (with a density of ~ 0.1 cm" 3 ). In addition, in the former case denser material (p 
~ 0.1 cm -3 ) from the progenitor star of the SN in about the first parsec of propagation is 
expected. 



4.3.2 Cannonball propagation in interstellar medium 

When a mass of magnetised plasma (plasmoid) moves through the ionised interstellar 
medium two limiting cases can be distinguished. 

1. If the plasmoid is impenetrable to the incoming charged particles (e.g. because of a 
homogenous, strong B field), so that all are effectively reflected, a collisionless bow-shock 
forms and the particles flow around the plasmoid without any sweep up of ambient matter. 

2. If the incoming charged particles can freely enter the plasmoid (e.g. because of a very 
inhomogeneous B field) no shock forms and all incoming particles are swept up. 
Because cannonball B-field is expected to be strong but very turbulent I think an interme- 
diate situation is most likely. The moving cannonball drives an ultrarelativistic collisionless 
shock into the interstellar medium, but most of the incoming particles (including the ones 
accelerated at this shock) are eventually swept up by the following cannonball. 

4.3.3 Numerical calculation of cannonball motion 

The following relativistic equations of motion for plasmoids ejected with an initial Lorentz 
factor of r ~ 300 and an initial energy of E ~ 10 53 ergs were numerically integrated using 
the following expressions fll3| : 

(10) dY/dm = -(r 2 - 1)/M 

(11) dm/dr = irR 2 p(r) 

(12) dM/dr = Tdm/dr 

m is swept up mass, p the ambient density and r the distance travelled by the cannonball. 
M is the cannonball mass, and was set to E/r initially. The cannonball radius R was 
calculated assuming that the plasmoid expands with a speed of v = 0.2 c/T in the observer 
frame until it reaches a radius of Rfi na i ~ 0.3 pc. No further expansion was assumed to take 
place afterwards, v is similar to the initial expansion in the plasma clouds ejected by GRS 
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1915+105 and also the speed of sound in a relativistic gas. 



As discussed in section 4,3.2 the plasmoid is assumed to sweep up all ambient matter, i.e. it 
grows in mass. The calculation does not include a treatment of the internal dynamics of the 
plasmoid, the matter is assumed to move with the plasmoids Lorentz factor at the time of 
sweep up. Therefore the presented results are only indicative and unrealistic, especially in 
the final stages of propagation when most of the initial kinetic energy has been transformed 
into turbulent plasma motion, B fields and CRs. 



For the ambient matter the two cases discussed in section 4.3.1 were assumed. In case 1 
the cannonball-plasmoid first moves in matter of 0.1 cm~ 3 for 1 pc (stellar wind), 1 kpc in 
10~ 3 cm~ 3 (super bubble) and then in 2 • 10 -4 cm -4 (halo). In case 2 the plasmoid moves 
in matter with a density of 0.1 cm -3 throughout. 

Some relevant parameters as a function of plasmoid travel distance are given in figs.^ and 
Case 1 (~ 85 % of all cannonballs, see section 4-3-j ) 



The plasmoid is slowed down non-ultra-relativistic speeds (r < 10) in the first ca. 40 pc. 
It then travels for a time t p (some hundreds of thousands years) with subrelativistic speed 
{(3 ~ 0.2). The number of cosmic-ray producing plasmoids at any time is given as t p /tj 
where tj is the mean interval of plasmoid producing events (estimated as 50 - 1600 years in 



section 4.1.2). I thus expect some hundreds to thousands of plasmoids in the halo that are 

still actively accelerating and releasing CRs. 

Case 2 (m 15 % of all cannonballs, see section 4-3. \ ) 



The cannonball Lorentzfactor V falls below 10 at about 10 pc and its speed (3 below 0.1 at 80 
pc. The cannonball does not travel far from its point of origin in a SN, i.e. it remains close 
to the Galactic centre in general. Consequently also all CRs produced by the cannonball 
will be released there, within the solar circle. 

These results were found to be very insensitive on the assumed initial value of F. However, 
if the initial energy E is assumed to be significantly smaller or the final confined plasmoid 
radius Rii m it significantly larger, all cannonballs range out before reaching the Galactic 
halo. 



5 Acceleration of cosmic rays by and in cannonballs 



5.1 lst-order Fermi acceleration up to the knee: ultrarelativistic shocks 



Heavens & Drury[34| - in a pioneering paper on acceleration in relativistic shocks - made the 
conjecture that the power-law index of particles accelerated in such shocks asymptotically 
approaches a universal index of 5 = —2.2 (eq.(Q)) with rising Lorentz factor. This has 
recently been confirmed by several researchers p9| , [|, 30 1. As this is exactly the source 



index of the "universal CR sources" in section 2.2, I propose that the bulk of Galactic 
hadronic CRs is accelerated at the ultra-relativistic shocks driven by cannonballs into the 
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Figure 6: The evolution of the Lorentz factor of a plasmoid with the properties described 
in the text as function of distance travelled in the observer frame. The full curve is for 
motion in a tenuous super-bubble medium (case 1), the dashed one for the normal warm 
interstellar medium (case 2). The dotted line shows the amount swept-up matter in the 
latter case in units of solar masses for case 2. 



interstellar medium. Extragalactic jets, which might be qualitatively similar objects, have 

been identified as prolific CR accelerators JlG|] , in particular their hot spots p7|. 

The accelerated particles are immediately swept up into the plasmoid (section |4.3.2 ) and 



remain confined there until it has slowed down to subrelativisitic speeds in the Galactic 
halo (section |5.2| ). In this scenario a spectrum extending down to small energies with the 
correct index is expected. The release mainly outside the solar circle leads to a spatial 
"source distribution" in agreement with the one expected for "universal sources" (section 
P). 

The maximal energy that can be reached in this acceleration process is limited by the 
extension of the shock and the interstellar magnetic field (the shock propagates in this 
medium) and is given aspCfl: 

(13) E max = 3 • W 15 eV{R/0.02pc){ZB/fiG){T/50) 

R is radius of the plasmoid when it has a Lorentz factor of T = 50. The preferred value 
of R=0.02 pc was chosen based on results of the numerical calculation discussed in section 
4.3.3] , Z is the charge of an accelerated particle. A precise calculation of this maximum 



energy would require a more detailed calculation of the plasmoid dynamics than I have 
performed, but the inserted numbers are correct to the order of magnitude. I identify E max 
with the "knee energy" Ek nee (sect |2.2| ). The observed knee seems to be a rather "sharp" 
feature in the primary CR spectrumQ. In my scenario there are ~ 1000 current sources 
for these CRs. It seems highly unlikely that they would be so uniform in their properties 
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Figure 7: The /3=v/c factor of a plasmoid propagating in a super-bubble environment 
(case 1) as a function of transversed distance in the observer frame. The dashed curve 
gives the elapsed time in the observer frame. 

that a sharp knee would be a natural feature]]. Therefore the sharpness of the knee remains 
unexplained. 

5.2 2nd-order Fermi acceleration above the knee: the Syrovatskii mech- 
anism 

While the plasmoids are "coasting" in the super-bubble (case 1 of section fl.3.3|) or normal 
interstellar medium (case 2) for a time t coas t ~ several 10 5 years (case 1) resp. hundreds 
years (case 2) with "subrelativistic" speeds internal turbulent plasma motions will further 
accelerate the hadronic "swept up" CR pool via the second-order Fermi mechanism. The 
coasting time scales t coas t for case 1 and 2 were estimated from the numerical calculation 
explained section 4.3.3 (see also fig.[|). 



Without turbulent plasma motion the swept-up CRs remain very efficiently confined in the 
plasmoid during the coasting phase except at the highest, near-knee energies. The diffusive- 
escape time scale from a sphere with radius R and a mean free scattering path A is given 



as trf=1.5 R 2 /cA [36]. Here A is a particles mean free path for scattering from magnetic 
inhomogeneities. I set A to its Larmor radius in this manuscript always (i.e. within the 
cannonball I always assume the Bohm limit). One then gets for the escape time from the 

7 I thank an anonymous referee for pointing this out to me 
8 defined here as 0.2 < (3 < 0.5) 
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plasmoid with the "confined" radius Rn m it and internal magnetic field strength B (section 



t d = 10 5 (i?i lmit /0.3pc) 2 x 
(14) (B/0.6G)/(E t /3 • 10 1E W)years « i coast 

At energies below Eknee td exceeds t coa st by a large margin. The rate of diffusive release - 



which determines the total CR-luminosity (section 4.1,3 ) - is then determined by the tur- 
bulent plasma motion in the plasmoid, and I do not try to determine it here. 
It is plausible that equipartition between turbulent motion, magnetic fields and cosmic rays 
holds in the confined plasmoid. The major energy loss then occurs via loss of cosmic-rays 
due to diffusion or release in bunches due to disturbances of the magnetic field at the plas- 



moid boundary. Syrovatskii|64f| has shown that in this situation equipartition is restored 
through particle acceleration and the released CR particles are distributed in energy ac- 
cording to a power-law spectrum with a universal differential index. In particular he shows: 
At energies that were not previously populated in the initial pool of CRs - i.e. in our case 
above the "knee" energy - the universal differential power-law index of the released CRs 
according to eq.(|3|) is 62 = —2.5. 

Moreover, it can be easily shown with the expressions in the appendix of Syrovatskii's paper 



At energies which were populated in the initial pool of CRs with an index of 5\=— 2.2 - i.e. 
in our case below the "knee energy" - the universal differential index of the released CRs 
is (again) —2.2. This "invariance" only holds for the index of —2.2 and is thus apparently 
fortuitous. 

Syrovatskii's results are independent of any details of the Fermi acceleration process 
because they follow from purely thermodynamic considerations. 

I thus obtain the source spectrum that we had derived as most likely for the "univer- 
sal source" class in section [2.2| , a power law with an index of —2.2 below the knee and 
—2.5 above the knee up to the highest observed energies / therefore propose that hadronic 
CRs mainly stem from cannonballs in the Galactic halo that release previously accelerated 
particles. 



5.2.1 Possible visibility of Galactic cannonballs in the radio range 

The expected radio luminosity of cannonballs in our Galaxy is rather smaller than one might 
expect at first sight. In their early ultrarelativistic stage relativistic beaming prevents the 
visibility of Galactic cannonballs. The emitted radiation is then visible only in external 
distant galaxies as a GRB afterglow. Galactic cannonballs become in principle visible only 
in the "coasting phase" , when no important beaming takes place and the cannonball accel- 
erates CRs via the Syrovatskii mechanism in the present scenario. In this phase nearly all of 
the CR acceleration has already taken place in the present scenario. Only a small fraction 
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/2nd-order (see below) of the total CR energy produced by a cannonball is accelerated and 
can contribute to the radio luminosity. 

A fundamental uncertainty in predicting the cannonball radio luminosity is the estimation 
of a factor f e / p , giving the CR energy injected into electrons relative to the one injected 
into nuclei in the coasting phase. Because of the high magnetic field in the cannonball, 
electrons have a very short lifetime in the plasmoid due to synchrotron losses and need to 
be injected during the "coasting phase" . It is a fundamental assumption of the Syrovatskii 
mechanism that the injection of new particles from the thermal pool has ceased, so that 



mainly already-relativistic particles are accelerated [64]. In this limit f e / p =0 and the radio- 
luminosity of cannonballs is small. However, it seems likely that this assumption holds only 
approximately and that f e / p is smaller than 1 by some not extremely large factor. The radio 
luminosity of a cannonball in the Galactic halo at 1.4 Ghz </>i.4GHz assuming a differential 
spectral index of the accelrated particles of -2.5 can be parameterized as follows: 

0i.4GHz ~ 3(£ tot /10 53 erg) x 
(-Rprog/lOOOyears) x 

(/2nd-ordcr/3 ■ 10~ 5 ) X 

(icoast/3 • 10 5 years) x 
(15) (20kpc/dcB) 2 /e/pJansky 

Here /2nd-order is the ratio of the energy injected into CRs in the second-order process 
relative to the total. The estimated value corresponds to the energy content of CRs above 
the knee relative to the total because in the present scenario CRs are reaccelerated to 
energies exceeding the "knee" energy in the "coasting phase" of the cannonballs. dcB is 
the distance of the cannonball. 

With a factor f e / p somewhat below 1 the expected radio brightness of cannonballs is in the 
sub-Jansky regime. A systematic search for nonthermal radio-sources with a high proper 
motion (on the order of 1 arcsecond/year) is a crucial test of the present scenario. There 
are 771000 radio sources the recent FIRST catalogue Q of which only 19 % have optical 
counterparts. No systematic search for high proper motions was performed in this data base. 



If one embarks on one, there are many possible technical reasons for false positives [ 72 ] . Even 
at a flux level of one Jansky not all radio sources in the FIRST catalogue are identified. 
Radio sources with flux levels near 1 Jy are frequently very faint high-z radio galaxies. 
In the absence of a proper- motion survey in the radio and/or systematic identifications of 
all radio sources with a limiting flux level below about a Jansky the absence of "obvious" 
radio counterparts to cannonballs cannot yet be an argument to reject the present scenario 
"a limite" . 

When the cannonball finally stops, it enters Sedov- Taylor phase. In this phase adiabatic 
losses probably prevent an important release of CRs and its brightness is difficult to predict. 
For an earlier very preliminary discussion about expected 7-ray fluxes of cannonballs see 
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5.2.2 Can second-order acceleration in the plasmoids reach ultra-high ener- 
gies? 

To reach the highest observed CR energies UH energies (i.e. energy of protons up at least 
Ej ~ 10 20 eV) in second-order Fermi acceleration three conditions have to be fulfilled [j36(] . 

1. The lifetime of the plasmoid "i C oast" (see beginning of section |5.2| ) has to exceed the 
time "t n " necessary for acceleration from the initial Ej to the final energy E^: 

(16) t C oast > hundreds of years > t n 

An analogous condition limits the maximum energy that can be reached in CR acceleration 
in SNR to about 100 TeV H. t n is given as ln(Ej/Ej) x t a . The second-order Fermi 
acceleration time scale is t a = A/2c/3 2 [36]. For a conservative estimate of t n I choose the 
A for the final maximum energy Ej = 10 20 eV. 

tn = 400(/3/0.5) 2 x 
ln{(E f /W 20 eV)/{Ei/3 ■ W 15 eV)) x 

(17) (Ef/10 20 eV)/(B/Q.6G)jeaxs 

Taking into account that this estimate is conservative I conclude that the condition ( |l6| ) 
can be surely met for protons if (3 > 0.5. 

2. The diffusive-escape time scale t^=l .5 R 2 /cA must be larger than the second-order 
acceleration time scale t a =A/2c/3 2 up to the highest energies, otherwise the particles escape 
before being accelerated. I obtain for protons: 

t d /t a ~ (/3/0.2) 2 ( J R limit /0.3pc) 2 x 

(18) (B/0.6G) 2 /(E/10 20 eV) 2 

The condition is thus fulfilled for /3 > 0.2. 

3. The synchrotron-loss time scale t s ~ k E B _ with k= 2e 4^ n ^ a ^ , has to remain larger 
than the acceleration time scale t a , otherwise acceleration is shut off. a is the pitch angle 
and is conservatively set to ir/2 below. 

t s = 3.9/((E/W 20 eV) ■ (5/0.6G) 2 )years w t a 

(19) m 3.1(E/10 20 eV)/((B/0.6G) x (/3/0.3) 2 )years 

The condition is seen to be fulfilled for (3 > 0.3. 

I conclude that there is no basic argument forbidding acceleration of UHE CRs in plasmoids. 
However, it seems difficult to accelerate protons to energies far above 10 20 eV with this 
mechanism. 
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6 Cannonballs as universal CR sources and observational 
facts 



6.1 UHE CRs from plasmoids in the Galactic halo: Greisen cutoff and 
anisotropics on a small angular scale 

If CRs originate in the Galactic halo, the absence of the "Greisen-Zatsepin cutoff" - that 
is otherwise made quite puzzling by the lack of extragalactic source candidates near the 
direction of individual UHE-CRs [65] - is readily understood. At a distance of typically 20 



kpc no significant absorption of UHE protons in the 3 K° background radiation is expected. 
There is some evidence for a small-scale clustering in the arrival directions of UHE CRs 



1 67]. The number of discrete sources necessary to explain the observed degree of clustering 



was estimated as ~ 400 [24] in agreement with the order of magnitude number of cosmic-ray 



producing plasmoids estimated in section 4.3.3. Consequently it is proposed to identify the 



subrelativisitic cannonballs in the Galactic halo as the sources responsible for the small-scale 



clustering [19 ] . 



6.2 Expected anisotropy of CR arrival directions from coasting cannon- 
balls 

6.2.1 At energies below the ankle 

Core-collapse supernovae are well known to occur predominantly within the solar circle ||. 
Therefore those cannonballs that range out in the dense disk-interstellar medium (corre- 
sponding to case 2 in section |4.3.3|) release all CRs at Galacto-centric distances smaller than 
the solar one. These CRs (estimated to be about 15 % of the total in section 4.3,3j ) will 



contribute to the anisotropy of CRs at all energies. We had seen in section 3.2 that about 
this fraction / (Eq.^) of CRs is required to be released by the "universal" sources within 
the solar circle. The cannonballs therefore remarkably "fit the bill" derived for "universal 



sources" in section 3.2 



6.2.2 At energies above the ankle 

The absence of significant large-scale anisotropies of arrival directions at UHE energies (sec- 
tion 2.3.2j ) is understood because an observer inside an emitting sphere sees no anisotropies. 



However about 15 % of all cannonballs are expected to release CRs at Galactic locations 
similar to the ones of SNe (section [4. 3. 3 ). Thus it is predicted that about 15 % of all CRs 



that propagate roughly on straight lines in the Galactic magnetic field (i.e. with energies 
well above the ankle) come from the general directions of the Galactic centre and disk. 
Presently the world-data sample at these energies is not large enough to confirm or rule 
out this prediction. 
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6.3 Galacto-centric distribution of the CR density 



I had discussed in section [3JJ that the very small observed gradient of CR density as a 
function of Galacto-centric radius strongly points towards CR sources mainly outside the 
solar circle. The cannonballs are expected to release ~ 15 % of all CRs within the solar circle 
(section 4.3.3; ), so a small non- vanishing gradient is expected. The data are in qualitative 
agreement with this expected "cannonball density distribution" but the presence of large 
systematic errors seem to make a further detailed test of the scenario in this area difficult. 



6.4 Predicted chemical composition of Cosmic Rays at all energies 

One can now calculate the chemical composition of cosmic rays. At low energies, below 
the knee, CRs are accelerated from the local interstellar medium and are thus expected to 
have abundances at a given energy per nucleus similar to the solar ones, as observed. I 
take the observed abundance at low energies as the starting point. In both first and second 
order Fermi acceleration identical power-law indices for all nuclei are expected. The knee 
(eq. (|l3|)) and ankle positions in energy are proportional to Z. With these assumptions the 



spectrum for all indices can be predicted (fig^). We saw in section 5.2.2 that acceleration of 
protons beyond x 10 20 eV is difficult in cannonballs and therefore assume an exponential 
cutoff at an energy/nucleus above Z • 10 20 eV. It is seen that in the universal scenario the 
chemical composition gets heavier near the knee - similar to what is expected in many other 
models. Near the ankle the composition becomes very light, similar to the expectations 
in some extragalactic scenarios for the origin of UHE CRs but finally, beyond 10 20 eV the 
composition becomes heavy again. 



7 Discussion 



7.1 Assumptions 

The major assumptions of the present model for CR origin are: 

1. The basic features of the "cannonball model" for GRBs [^0| are correct. 

2. The CR propagate diffusively in a relatively large Galactic halo (radius ~ 20 kpc). The 
energy dependence of the diffusion coefficient is Kraichnanian at all energies. 

3. The total energy of the ejected cannonballs is ~ 10 53 ergs, about an order of magnitude 
higher than often assumed. 

4. After an initial expansion the cannonball remains confined to a radius of about 0.3 pc, 
due to ram-pressure and magnetic effects until it has slowed down to speeds below about 
0.1 c. 

While there are theoretical and observational arguments in favour of assumption 1 and 2, 
assumptions 3 and 4 have no better motivation presently than to make the model work. 
Additional minor assumptions are discussed in the text. 
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Energy In eV 



Figure 8: Upper panel. The intensity of 5 chemical groups (in the order of falling intensity 
at low energies: p (dashed), He (dotted), CNO group (dash-dot), Ne-S group (dashed), Cl-Fe 
group(full)) in m -2 s _1 sr _1 TeV 1 ' 7 as a function of energy /nucleus. The uppermost line 
is the sum of all components in the same units. The composition for low energies below 
the knee is taken from Wiebel[74|. Above the knee the prediction of the present model 
is shown. Lower panel: the ratio of the light element intensity (p + He) to the total 
intensity, as derived from the upper panel. 



7.2 Achievements 

The location of CR sources at all energies mainly in the Galactic halo and with a fraction / 
~ 15 % near the Galactic centre expected in our model allows to understand the following 
observations more naturally than many other models: 

1. The absence of the Greisen-Kuzmin cutoff together with the absence of a large scale 
anisotropy at energies above about 4 x 10 19 eV and extragalactic counterparts in the di- 
rection of the highest-energy CRs. 

2. The observed anisotropy of CR with energies between about 100 TeV and 10 18 eV quan- 
titatively. 

3. The Galacto-centric distribution of low-energy hadronic CRs inferred from 7-ray astron- 
omy. 

The origin of hadronic CRs in about 1000 cannonballs speeding through the interstellar 
and interhalo medium allows to understand: 

4. The possible indication for small-scale clustering is understood as being due to emission 
from individual cannonballs. 
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5. The origin of the "knee feature" in the hadronic CR spectrum and the form of the 
source-spectrum and higher and lower energies. The observed "sharpness" of the knee re- 
mains unexplained, however. 

6. Finally the idea of a "universal-source class", accelerating CRs at all energies, leads 
to a very natural understanding of the energy for the "ankle" feature in the spectrum of 
Galactic CRs. 



7.3 Predictions 

The most characteristic prediction of a universal origin of hadronic CRs at all energies 
is the variation of chemical composition above the knee shown in fig.(||). Evidence for 
a Kolmogoroff rather than Kraichnan energy dependence of the diffusion coefficient for 
cosmic rays (e.g. from the ongoing experimental studies on reacceleration) would rule out 
the "universal scenario" presented here. If "cannonballs" accelerate all CRs above the 
ankle, on the order of 1000 spatially discrete sources of UHE CRs must be found by the 
next generation of UHE-CR detectors. About 15 % of these clusters must be preferentially 
spatially distributed like core-collapse SNe (i.e. near the Galactic disk and centre). 
The counterparts of these sources must be very compact plasma clouds with subrelativistic 
proper-motions. The present model stands or falls with the existence of counterparts to 
the cannonballs in various regions of the electromagnetic spectrum. 

Acknowledgements I thank Maria Diaz for reading the manuscript and Arnon Dar 
and Alvaro De Rujula for many helpful explanations of the cannonball model. I am grateful 
to an anonymous referee for helpful criticism. The author is supported by a Heisenberg 
fellowship of the DFG. 



References 

[1] Armstrong J.W., Cordes J.M., Rickett B.J., 1981. Density power spectrum in the 
local interstellar medium, Nature 291,561-564. 

[2] Arqueros et al. (HEGRA collaboration), 2000. Energy spectrum and chemical com- 
position of cosmic rays between 0.3 and 10 PeV determined from the Cherenkov-light 
and charged particle distribution in air showers, A & A 359,682-694. 



[3] Becker R.H., et al., 2001. FIRST, |http://sundog.stsci.edu/ 



[4] Bednarz J., Ostrowski M., 1998. Energy Spectra of Cosmic Rays Accelerated at 
Ultrarelativistic Shock Waves, PhRvL 80,3911-3914. 

[5] Berezinskii V.S., Bulanov S.V., Dogiel V.A., Ginzburg V.L., Ptuskin V.S. , 1990. 
Astrophysics of cosmic rays. North Holland, Amsterdam. 



26 



[6] Biermann L., Davis L., 1958. On the origin of Cosmic Rays during the early part of 
the Evolution of our Galaxy, Z.Naturforsch. 13a, 909-915. 



[7] 
[8] 

[9 

[10 
[11 



[12; 
[13; 

[14 



[15 



[16 



[17 



[is; 



[19 



[20 



[21 



Biermann P.L., 1995. The origin of cosmic rays, PhRvD 51,3450. 

Bloemen H., 1993. Diffuse Galactic Gamma- Ray Emission Above 1 MeV, AIP Conf. 
Proc. 278,409. 

Burbidge G., 1955. Halo of radio emission and the Origin of Cosmic Rays, Phys.Rev. 
101,906-907. 

Burbidge G., 1962. The origin of Cosmic Rays, Prog. Theoret. Phys. 27,999. 



Cen R.,1998, Supernovae, pulsars and gamma-ray bursts: a unified picture, astro- 



ph/9809022| , ApJL. 



Cesarsky C.J. ,1980. Cosmic Ray propagation in the Galaxy, ARA&A, 18, 289-319. 

Chiang J., Dermer C.F., 1998. Synchrotron and SSC Emission and the Blast-Wave 
Model of Gamma-Ray Bursts, astro-ph /980333S . 

Clay R.W., Donough M.-A., Smith A.G.K., Dawson B.R., 1998. Anisotropies 
and the power requirements for Galactic Cosmic Rays, PASA 15, available at 



tittp : / / www . atnf . csiro . au/pasa/ 152/ 



Compton A.H., Getting I.A., 1935. An apparent Effect of Galactic Rotation on the 
Intensity of Cosmic Rays, Phys.Rev. 47,817-821. 

Connell J.J., Simpson J. A., 1999. Testing the role of Cosmic Ray reacceleration, Proc. 
26th ICRC,vol.3,p.33-26. 



Dar A., 1998. Cosmic rays and gamma-ray bursts from microblazars, astro- 



ph/9809163 



Dar A., 1999. Are extragalactic Gamma Ray Bursts the Source of the Highest Energy 



Cosmic Rays?, astro-ph/9901005. 

Dar A., Plaga R., 1999. Galactic 7-ray bursters - an alternative source of cosmic rays 
at all energies, A&A 349,259. 

Dar A., De Riijula A., 2000a. A cannonball model of gamma-ray bursts: superluminal 
signatures, astro-ph / 0008474 , A & A in the press. 

Dar A., De Riijula A., 2000b. The Cannonball Model of GRBs: Temporal and Spec- 
tral Properties of The Gamma Rays, astro-ph/0012227, submitted A Sz A. 



27 



[22] Dermer O, 2000. Neutrino, Neutron, and Cosmic Ray Production in the External 
Shock Model of Gamma Ray Bursts, |astro-ph/0005440 , published sparately in: ApJ 
574 in the press; and ApJ 556, 479. 

[23] De Young D.S., Axford W.I., 1967. Inertial confinement of Extended Radio Sources, 
Nature 216,129. 

[24] Dubovsky S.L., Tinyakov P.G., Tkachev I.I. ,2000. Statistics of clustering of ultra-high 
energy cosmic rays and the number of their sources, astro-ph/0001317, Phys. Rev. Lett. 
85, 1154. 

[25] Erlykin A. D., Wolfendale A.W., Zhang L., Zielinska M., 1996. The gradient of cosmic 
ray protons in the outer Galaxy, A&ASS 120,397. 

[26] Erlykin A. D., Wolfendale A.W., 2000. The origin of PeV cosmic rays, A&A 356,L63- 
L65. 

[27] Fenimore E.E., Ramirez-Ruiz E., 2000. Redshifts for 220 BATSE Gamma-Ray Bursts 
Determined by Variability and the Cosmological Consequences, |astro-ph/0004176| . 

[28] Fishman, G. J. & Meegan, C. A. A. 1995. Gamma-Ray Bursts, ARA&A 33,415. 

[29] Fryer C.L., Woosley S.E, 1999. Formation Rates of Black Hole Accretion Disk 
Gamma- Ray Bursts, astro-ph/9904122 . 

[30] Gallant Y.A., Achterberg, A. ,1999. Ultra-high-energy cosmic ray acceleration by rel- 
ativists blast waves Jastro-ph/9812316l v2, MNRAS 305, L6. 

[31] Ginzburg V.L., 1957. The origin of cosmic rays, Prog. Elem. Particle Cosm. Ray 
Phys. 4,339 

[32] Ginzburg V.L., Ptuskin V.S., 1976. On the origin of cosmic rays: Some problems in 
high-energy astrophysics, RvMP 48,161-189. 

[33] Hayashida N. et al. (AGASA coll.), 1998. The Anisotropy of Cosmic Ray Arrival 
directions around 10 18 eV, |astro-ph/9807045| , APh 10, 303-311. 

[34] Heavens A.F., Drury, L.OC.,1988. Relativistic shocks and particle acceleration, MN- 
RAS 235,997-1009. 

[35] Heinz S., 2000. Beams, Bursts, Bubbles and Bullets: Relativistic Out- 
flows in Astrophysics, PhD thesis University of Colorado, chapter 5.6.3, 
http : // jilawww .color ado.edu/ www / sro/ thesis /heinz/ index. htm\ 



[36] Hillas A.M., The origin of Ultra-High-Energy Cosmic Rays, 1984. ARA & A, 22,425- 
444. 



28 



[37] Hunter S.D et al., 1997. EGRET Observations of the Diffuse Gamma- Ray Emission 
from the Galactic Plane, ApJ 481,205. 

[38] Kembhavi A.K., Narlikar J.V.,1999. Quasars and active galactic nuclei. Cambridge 
University Press, Cambridge. 

[39] Kirk J.G., Guthmann A.W., Gallant .A., Achterberg, A.. 2000. Particle acceleration 
at ultrarelativistic shock: a eigenfunction method, astro-ph/0005222 , ApJ 542,235. 

[40] Lamb, D.Q., 2000. Implications of recent observational dicoveries for the nature and 
origin of gamma-ray bursts, |astro-ph/0005028| . 

[41] Lanzetta K.M., Bowen D.V, Tytler D.,Webb J.K., 1995. The Gaseous Extent of 
Galaxies, ApJ 442,538-568. 

[42] Lilly, S.J. et al. 1996. The Canada-France Redshift Survey: The Luminosity Density 
and Star Formation History of the Universe to Z approximately 1, ApJ 460, LI. 

[43] Lingenfelter R.E., Higdon J.C, Ramaty R., 2000. Cosmic Ray Acceleration In Super- 
bubbles ad the Composition of Cosmic Rays, astro-ph/0004166 . 

[44] Linsley J., 1983. Spectra, anisotropies and composition of cosmic rays above 1000 
GeV, Proc. 18th ICRC Bangalore 12,135. 

[45] Longair M.S. ,1997. The high-energy radiation of active galactic nuclei, In: G. Munch 
et al. (Ed.), The Universe at Large. Cambridge University Press, Cambridge, p. 218. 

[46] Longmire C.L.,1963. Elementary Plasma Physics, Interscience (New York, London) 
p.68-72. 

[47] Loveday, J. et al. 1992. The Stromlo-APM Redshift Survey. I - The luminosity func- 
tion and space density of galaxies, ApJ, 390, 338. 

[48] Milgrom M., Usov V., 1996. Gamma-ray bursters as sources of cosmic rays, Astropart. 
Phys. 4,365. 

[49] Morrison P., 1957. On the origins of Cosmic Rays, RvMP 29,235. 

[50] Nagano M., Watson A. A., 2000. Observations and implications of the ultrahigh- 
energy cosmic rays, RvMP 72,689-732. 



[51] Paczyhski B., 1997. Are gamma-ray bursts in star forming regions': 



astro- 



ph/9710086t ApJ 494,L45. 

[52] Paczynski B., 1999. Gamma-Ray Burst - Supernova Relation, |astro-ph/990904£ , In: 
M. Livio, N. Panagia and K. Sahu (eds.), Proc. of the Space Telescope Science In- 
stitute 1999 May Symposium: "Supernovae and Gamma Ray Bursts; The Largest 
Explosions Since the Big Bang". Cambridge University Press, Cambridge, p.l. 



29 



[53] Plaga R., 1998. An extragalactic "flux trapping" origin of the dominant part of 
hadronic cosmic rays? A & A 330,833. 



[54] Plaga R., 2001. The cepheid-like relationship between variability and luminosity ex- 
plained within the "cannonball model" of Gamma-Ray bursts, |astro-ph/001206C , 
A&A 370,351. 

[55] Plaga R.,de Jager O. & Dar A., 1999. Are Galactic Gamma- Ray Bursters the Main 
Source of Hadronic Non-Solar Cosmic Rays at all Energies? Proc. 26th ICRC, Salt 
Lake City, 4,353. 

[56] Pugliese G., Falcke H., Wang Y. & Biermann P.L., 2000. The jet-disk symbiosis 
model for Gamma Ray Bursts: cosmic ray and neutrino background contribution, 
astro-ph/0003025| , A&A 358,409. 

[57] Rachen J., Biermann P.L., 1993. Extragalactic ultra high energy cosmic rays, I. Con- 
tribution from hot spots in FR-II galaxies, A&A 272,161. 

[58] Ramirez-Ruiz E., Fenimore E.E., 1999. Pulse width evolution in gamma-ray bursts: 
evidence for internal shocks, astro-ph/9910273] , ApJ 539,712. 

[59] Scalo J., Wheeler J.C., 1999. Astrophysical and Biological Implications of Gamma 
Ray Burst Properties, astro-ph /99 1 2564| . 

[60] Schmele D., 1998. Suche nach raumlichen Anisotropien der kosmischen Strahlung im 
Energiebereich oberhalb von 20 TeV, PhD thesis, University of Hamburg, available 



on http: //www.desy . de~schmele/D0KT0R/ 

[61] Spitzer L., 1990. Theories of the hot interstellar gas, ARA&A 28,71-102. 

[62] Strong A.W., Mattox J.R., 1996. Gradient model analysis of EGRET diffuse Galactic 
gamma-ray emission, A&A 308, L21. 

[63] Strong A.W., Moskalenko I.V., 1998. Propagation of Cosmic-Ray Nucleons in the 
Galaxy, |astro-ph/980715q v2, ApJ 509,212. 

[64] Syrovatsky S.I., 1961. Spectrum of Galactic and Solar Cosmic Rays, Sov.Phys. JETP 
13,1257-1260. 

[65] Takeda M. et al. (AGASA coll.), 1998. Extension of the Cosmic-Ray Energy Spectrum 
Beyond the Predicted Greisen-Zatsepin-Kuzmin Cutoff, Phys. Rev. Lett. 81,1163; 



astro-ph/9807192 



[66] Tkaczyk W., 1995. The Cosmic Ray Spectrum and irregularities in the Galactic 
magnetic field, Proc. 24th ICRC, Rome, 3,293. 



30 



[67] Uchihori Y., Nagano M., Takeda M., Teshima M., Lloyd-Evans J. & Watson A. A., 
1999. Cluster Analysis of Extremely High Energy Cosmic Rays in the Northern Sky, 
astro-ph/9908T93| , APh 13,151-160. 

[68] Vietri M., 1995. On the acceleration of Ultra High Energy Cosmic Rays in Gamma 
Ray Bursts, ApJ 453,883. 

[69] Wang L., Howell D.A., Hoflich P., Wheeler J.C., 1999. Bi-polar Supernova 
Explosions ,|astro-ph/9912033 . 



[70] Waxman E., 1995. Cosmological Gamma Ray Burst and the highest energy Cosmic 
Rays, PhRvL 75,386. 

[71] Weiner B.J., Williams T.B., 1999. Detection of Ha Emission from the Magellanic 
Stream: evidence for an Extended Gaseous Galactic Halo ,AJ 111,1156-1163. 

[72] White R.L., 20001. Private communication, January 5. 

[73] Wibig T., Wolfendale A., 1999. The anisotropy of arrival directions of cosmic rays 
above W 17 eV, J.Phys.G 25,2001-2009. 

[74] Wiebel B., 1994. Chemical composition in high energy cosmic rays, Report WUB 
94-08, available on http:// wpos6.physik.uni-wuppertal.de: 8080 
/Public/papers-public .html 

[75] Woosley S.E.,1993. Gamma-ray bursts from stellar mass accretion disks around black 
holes, ApJ 405,273. 



31 



